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ABSTRACT 


Saudi Arabia is enriched by many geothermal resources that are located along the western coastal part of 
the Red Sea, in the form of a number of hot springs and many volcanic eruptions. Wadi Al-Lith is 
considered one of the most promising geothermal targets that encounter many hot springs with a good 
surface temperature upto 95 °C. This paper aims mainly to evaluate the geothermal potential of the main 
hot spring at Wadi Al-Lith area (Ain Al Harrah). The available remote sensing images are analyzed and a 
number of 2D electric resistivity profiles are interpreted to delineate the surface geological lineaments 
and the subsurface structural elements that control the movement of the thermal water. It is found that 
the main surface lineaments are structurally oriented along NNE-SSW and NE-SW directions with a 
frequency percentage of 52% and an average lineament length of 835 m. Furthermore the subsurface 
structural elements, as inferred from the interpretation of the geophysical 2D electric profiles, have 
assumed the same directions beside the NW-SE direction. 

The characteristics of thermal water of the hot spring are indicated through analyzing the major and 
minor elements of some collected water samples. Geo-thermometers are applied to estimate, subsurface 
temperature, heat flow and discharge enthalpy. These parameters are found to be 136 °C, 183 mW/M? and 
219 kJ/kg, respectively. An estimate of the geothermal reserve using the volumetric method, gave total 
stored heat energy of 1.713 x 10” J (rock and fluid) and a geothermal reserve potential of 26.99 MWt. 

It appears from our research that the estimated energy is quite enough for operating a medium scale 
power plant that utilizes low boiling point fluid (Kalina Cycle) for limited electricity production, beside 
other low-temperature applications (district heating, green houses and medical therapy). 

© 2013 Elsevier Ltd. Open access under CC BY-NC-SA license. 
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1. Introduction 


Geothermal power is among the renewable energy sources which 
come from natural resources such as sunlight, wind, rain, tides, and 
biomass which are renewable. About 16% of global final energy 
consumption comes from renewables, with 10% coming from tradi- 
tional biomass, which is mainly used for heating, and 3.4% from 
hydroelectricity. New renewables (small hydro, modern biomass, 
wind, solar, geothermal, and biofuels) accounted for another 3% and 
are growing very rapidly. The share of renewables in electricity 
generation is around 19%, with 16% of global electricity coming from 
hydroelectricity and 3% from new renewables [1,2]. 

The most critical factor for the classification of geothermal energy 
as a renewable energy scale is the rate of energy recharge. In the 
exploitation of natural geothermal systems, energy recharge takes 
place by advection of thermal water on the same time scale as 
production from the resource. This justifies our classification of 
geothermal energy as a renewable energy resource. In the case of 
hot, dry rocks, and some of the hot water aquifers in sedimentary 
basins, energy recharge is only by thermal conduction; due to the slow 
rate of the latter process, however, hot dry rocks and some sedimen- 
tary reservoirs should be considered as finite energy resources [3]. 

Geothermal power is cost effective, reliable, sustainable, and 
environmentally friendly, but has historically been limited to areas 
near tectonic plate boundaries. Recent technological advances 
have dramatically expanded the range and size of viable resources, 
especially for applications such as home heating, opening a 
potential for widespread exploitation. Geothermal wells release 
greenhouse gases trapped deep within the earth, but these 
emissions are much lower per energy unit than those of fossil 
fuels. The International Geothermal Association (IGA) has reported 
that 10,715 MW of the online geothermal power in 24 countries, 
are generating 67,246 GWh of electricity in 2010. This represented 
a 20% increase in geothermal power online capacity since 2005. 
IGA expects that this will grow to 18,500 MW by 2015, due to the 
large number of projects presently under consideration, often in 
areas previously assumed to have little exploitable resource [4,5]. 

In Saudi Arabia there is a big gap of knowledge and information 
concerning the potential of geothermal resources. A more detailed 
and much more comprehensive work is needed to study the 
geothermal resources in Saudi Arabia, especially those encountered 
in Southwestern parts, from one hand, and to evaluate their 
economic reserves for possible energy production and domestic 
application, from the other hand. Large number of exploration 
methods and technologies should be initiated in order to reach these 
objectives. Many of these methods are in current use and have 
already been widely experimented in other sectors of research [6]. 

The geothermal recourses of Saudi Arabia can be categorized as 
the following: 


1. Low enthalpy resources (sedimentary aquifers): represented by 
deep-seated aquifers encountered in thick sedimentary basins 
in the eastern part of the Kingdom. These resources are 


confined and geothermal potential is represented by the 
normal geothermal gradient. It can be accessed only by deep 
drilled wells (i.e. oil wells). 

2. Medium enthalpy resources (hot springs): encountered along the 
western and southwestern coastal parts and represented by the 
shallow hot springs of hot surface water. These resources are 
unconfined targets with direct accesses to the subsurface hot 
anomalies through an open network of active faults and 
fractures (structure control). 

3. High enthalpy resources (basaltic lavas, Harrats): Saudi Arabia 
has approximately 80,000 km? of lava fields, known as Harrats. 
They are represented by volcanic eruptions, mainly basaltic in 
composition, that extend along the coastal part of the Red Sea 
at the western of Saudi Arabia. Harrats of Khaybar and Rahat 
are believed to be the best in terms of high heat flow and 
enthalpy [6]. 


The first geothermal work done at Saudi Arabia was that of El 
Dayel, 1988 [7]. He focused mainly on studying the hydro-chemical 
characteristics of the hot springs at Jizan and Al-Lith areas. Some 
water samples from the hot springs were analyzed. Other inter- 
ested work regarding the petrography, petrochemical composition 
and mineralogy of the basaltic volcanic eruptions, was done by 
Roobol et al. and Pint et al. in the period from 1992 to 2007. A 
more detailed and recent work is done by Lashin and Al Arifi, 2012 
[8]. They focused on giving a preliminary investigation of the 
geothermal resources encountered in the southwestern parts of 
the Saudi Arabia (Jizan area), based on the interpretations of the 
Landsat images and chemical analyses. 

In general, a few studies have presented the characteristics and 
the potentiality of the geothermal resources and its utility for 
some locations along the western coastal area of Red Sea of Saudi 
Arabia [8-14]. This work presents an initial evaluation of the 
geothermal energy potential of Wadi Al-Lith area, and provides a 
basis for a preliminary assessment of geothermal resources in the 
whole region. 


2. Geologic setting 
2.1. Rock units 


Wadi Al-Lith is represented by a catchment area that is 
considered an integral part of the Arabian Shield. It extends 
from the western coast of the Red Sea to high Mountains in the 
east. The area is mainly covered by metavolcanic rocks, metasedi- 
ments and late Proterozoic plutonic rocks (Fig. 1). Four major rock 
units are found covering about 91% of Wadi Al-Lith catchment 
area. These rock units are summarized, from youngest to oldest, as 
follows: 


1- Quaternary, sand, gravel and silt: This unit represents 18% of the 
total catchment area of Wadi Al-Lith. Eolian sand-dune fields, 
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Fig. 1. Geologic map of Wadi Al-Lith area showing Ain Al Harrah hot spring. 


mobile sheet sand and silt deposits cover a large proportion of 
the coastal plain in the southern part of the study area. Barrier, 
longitudinal, and barchans dunes are present locally, and most 
of the area is covered by sand deposits. 

Late- to post-tectonic granitic rocks: These rocks are part of the 
Proterozoic intrusive rocks that cover most of the western part 
of the study area. They represent 13% of the total basin area of 
Wadi Al-Lith. The plutonic rocks range in composition from 


serpentinite to syenite, although diorite, tonalite, granodiorite, 
and monzogranite predominate. 


3- Lith suite, Khasrah complex, diorite and gabbro: These rocks 


contain mafic to intermediate plutonic rocks that crops out in a 
wide northeast trending belt in the western part of Wadi Al- 
Lith basin. It covers almost 37% from Wadi Al-Lith basin area. 


4- Baish and Baha groups: Many rocks assigned to the undivided 


Baish and Bahah groups crops out in the eastern part of Wadi 
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Al-Lith basin as “basalt and dacite”, and “biotite-hornblende The main structural elements which are prevailing along the 
schist, and amphibolite” covering almost 23% from Wadi Al-Lith western coastal part (Red Sea) of the Saudi Arabia, including the 
basin area. study area, can be classified as: 

Other rock units are found constituting about 9% of wadi Al- è Precambrian northwest-southeast to east-west compression 
Lith basin. The most important of these rocks are the metagabbro, that are reactivated by tensional stress field during the Tertiary 
gabbro dikes, syentectonic granitic rock, granodiorite, granite and evolution of the Red Sea rift. 
amphibolite, as shown on the geological map (Fig. 1). © Tertiary northeast-southwest tension. 
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These structural elements are believed to be the main channels 
and pathways through which the thermal water rises to the surface. 


2.2. Topography 


Fig. 2 exhibits the digital elevation model of Wadi Al-Lith basin, 
which shows a wide variation in the topography of the area. 
A maximum depth of 2600 m (above sea level) is recorded at 
the highly elevated mountainous area in the northern part of Wadi 
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Al-Lith basin, while an elevation of less than 1 m is recorded for 
the sea warded coastal plains in the southern part of the area. 

According to its surface elevation from the sea level, Wadi 
Al-Lith basin can be divided into three zones, as follows: 


1- A high mountainous upstream zone: the elevation ranges 
from 500 m to 2600 m above the sea level and is characterized 
by high degrees of slopes and rugged area. It includes a number 
of tributaries and sub-wadies such as; Wadi Birayn, Wadi 
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Fig. 3. Landsat RGB image of Wadi Al-Lith catchment area. 
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Salibah and Wadi Magqsa which join the mainstream of Wadi 
Al-Lith. 

2- An intermediate zone of small hills and mountains: the eleva- 
tion of this zone ranges from 100 m to 500 m above sea level. It 
represents the midstream of the Wadi Al-Lith and includes the 
main tributaries and wadies sloping from up-stream. It is 
characterized by moderate gradient and encounters the major 
waterways of the Wadi Al-Lith. 

3- The coastal flood plain down-stream zone: this zone contains 
alluvial deposits with an elevation range from 100 m to nearly 
zero at the shoreline of the Red Sea. 


3. Methodology 


The utilized methods in this work are mainly applied to detect 
the following issues: 


1- The surface and subsurface structural elements that control the 
movement of the ascending thermal water. 

2- The petro-thermal characteristics of thermal water. 

3- The possible geothermal energy potential and its reserves. 


The analysis of surface and subsurface structural elements are 
conducted through analyzing the available remote sensing images 
and conducting a number of geophysical 2D electric profiles. 


Table 1 
The chemical analyses of Ain Al Harrah hot spring. 


3.1. Landsat RGB images 


The remote sensing imagery is useful in mapping lithology and 
offers an efficient tool for detecting geological structures, such as 
lineaments. In the present study, extraction of lineaments is enhanced 
through analyzing the available Landsat RGB images with 30m 
resolutions (Fig. 3). PCI-Geomatica software is used as a powerful tool 
for extracting the surface geological structures in the study area. 


3.1.1. Gridding and mapping of lineation 

The lineation data that are extracted from the analysis of 
Landsat images are feed to specific software (RockWorks Utilities) 
to create a grid file that models all aspects of lineation frequency, 
cumulative lengths, and/or intersections, with a variety of weight- 
ing options. Three main maps were constructed: 


1. Lineament intersection map: It illustrates the intersection grid 
map in which each cell represents the number of intersections 
of the lineations that occur within the cell. 

2. Lineament length map: It is a length grid map, in which each cell 
represents the total added lengths of all of the lineations that 
pass through this cell (total length of each lineation is 
included). 

3. Lineament frequency map: It is a lineation frequency grid map 
that includes the number of lineations that lie in the cell. 


On the other hand, a 2D geophysical electric resistivity method, 
is applied to investigate the subsurface orientation of the 


Locality Wadi Al-Lith Sampling date 28 May 2011 TDS 2180 ppm 
Name & No. Ain Al Harrah Type of aquifer Hot spring PH 78 
Longitude 40°28/16.4! Temperature <90°C EC 3633 pS cm! 
Latitude 20°27/41!! Elevation 167.0 m 
Chemical composition (mg/l) 
SiO2 Na* Kt Catt Mg** S04 NH4 HCO;37 cd F 
72.69 403.5 25.2 260 0.461 225.4 * 18.481 845 2.53 
Chemical composition (107° mg/l) 
Li Al Zn Fe Cu Mn Ni Co Ba B 
90.5 6.9 3.39 647 3.66 1.35 3.8 0.393 80.9 378 
Chemical composition (ppb) 
Rb Cr V Mo Sr U Th Not detected 
68.9 87.4 28.9 34.3 4160 0.216 0.0217 Over range 

Table 2 

The chemical analyses of water well close to the hot spring. 
Locality Wadi Al-Lith Sampling date 28 May 2011 TDS 1127 ppm 
Name & No. 28 Type of aquifer Water well PH 75 
Longitude 20°17/32! Temperature 34°C EC 1879 uS cm! 
Latitude 40°28/36.1!/ Elevation 60.0 m 
Chemical composition (mg/l) 
SiO2 Nat K+ Catt Mg++ S077 NH, HCO3~ cr F 
39.28 212.6 9.34 103 37.8 340 38.872 185.9 0.61 
Chemical composition (107° mg/l) 
Li Al Zn Fe Cu Mn Ni Co Ba B 
1.68 9.21 2.93 270 1.91 0.638 4.06 0.294 27.4 505 
Chemical composition (ppb) 
Rb Cr Vv Mo Sr U Th Not detected 
0.587 56.8 41 23.2 1630 2.6 Over range 
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Fig. 5. Rose diagram showing the structure lineaments of (A) Wadi Al-Lith area and (B) Ain Al Harrah hot spring area. 
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Table 3 
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The directional trends of the extracted lineaments in the Wadi Al-Lith area. 


Lineaments structure trend 


Degrees range Frequencies A 


N-S 
NNE-SSW 
NE-SW 
ENE-WSW 
E-W 
ESE-WNW 
SE-NW 
SSE-NNW 


(0-11.25), (168.75-191.25) and (348.75-360) 1644 
(11.25-33.75) and (191.25-213.75) 2297 
(33.75-56.25) and (213.75-236.25) 2230 
(56.25-78.75) and (236.25-258.75) 1683 
(78.75-101.25) and (258.75-281.25) 839 
(101.25-123.75) and (281.25-303.75) 329 
(123.75-146.25) and (303.75-326.25) 339 
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(146.25-168.75) and (326.25-348.75) 815 ; 


Table 4 


The directional trends of the extracted lineaments of Ain Al Harrah area. 


Lineaments structure trend 
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Fig. 6. Lineation intersection grid map of Wadi Al-Lith basin including Ain Al Harrah area. 
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Fig. 7. Lineation length grid map of Wadi Al-Lith basin including Ain Al Harrah area. 


geothermal targets and detecting their possible lateral extension 
as well as the structural elements which may control the subsur- 
face movement of the arising geothermal water. 


3.2. Geophysical surveys 


The geothermal potential can be explored and investigated 
using a variety of geophysical and geochemical (geo-thermo- 
meters) techniques which are normally used in geo-technical 
and geological investigation as well as in the exploration of 
hydrocarbons in the oil industry [14-22]. Such techniques provide 
the necessary information for constructing a realistic model of the 
geothermal system and assessing the potential of the resource. The 
electric resistivity is related to various geological parameters such 
as the mineral and fluid content, porosity and degree of water 
saturation in the rock [23]. 

A number of 2D electric profiles are performed at the location 
of the hot springs in Wadi Al-Lith area. The Syscal-R1 system (IRIS 
instruments) with 72 multi-electrodes system, with 4 x 18 elec- 
trode cables is used for conducting the field surveys using different 
electrode spacing (1m, 3m, 5m and 10m) arranged along a 


profile. Wenner-Schlumberger arrangement is used and the sur- 
veyed sections are taken very close to and crossing the hot spring. 
RES2DINV software is used for processing, correcting and inter- 
pretation of electric resistivity data. It can operate large datasets 
collected by a large number of electrodes and can account for the 
topographical effect along the survey lines. Finally, it generates 
two-dimensional resistivity model for subsurface from electrical 
imaging surveys [23]. 


3.3. Petro-thermal characteristics of thermal water 


The petro-thermal characteristics of thermal water of Ain Al 
Harrah hot spring at Wadi Al-Lith area are mainly enhanced through 
analyzing water samples collected from the hot spring and a very 
close water well. These water samples are analyzed using the 
ICP-MS spectrometer. The analysis was performed by Perkin-Elmer 
Sciex instruments multi-element ICP-MS spectrometer (type 
ELAN9000) equipped with standard torch, cross flow nebulizer 
and Ni sampler and skimmer cones. Tables 1 and 2 show the 
chemical analyses (major, minor and trace elements) of two water 
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Fig. 8. Lineation frequency grid map of Wadi Al-Lith basin including Ain Al Harrah area. 


samples taken from Ain Al Harrah hot spring and the neighboring 
water well. 


3.3.1. Geo-thermometers 

Geo-thermometer is a chemical substance or isotope in water 
that may be used to predict the subsurface temperature in 
geothermal systems. The basic assumption behind using the geo- 
thermometers is the chemical equilibrium between the substance 
(isotope) and mineral (s) in the reservoir. Theoretically, any cation 
ratio and any aqueous species concentration can be used as a geo- 
thermometer as long as equilibrium prevails. A temperature 
equation for a geo-thermometer is a temperature equation for a 
specific equilibrium constant referring to a specific mineral- 
solution reaction [24-29]. 

A selected number of geo-thermometers (Quartz, Chalcedony, 
Na, K, Ca geo-thermometers and combinations of them; i.e. Na-Ca 
and Na-K-Ca) are applied for determining the subsurface forma- 
tion temperature, discharge enthalpy and heat flow. These geo- 
thermometers have a low to medium temperature range of 25- 
250 °C that fit with the thermal properties of the hot spring under 


study. Eqs. (1)-(7) demonstrate the mathematical expressions of 
these geo-thermometers [24-29]: 

(a) Silica geo-thermometers 

Quartz 


TCO = 579302 g;-273.15 (Range of °C : 25-250) (1) 


T2(°C) = 42.198 + 0.288310z-3.6686 x 10-4Qz? + 3.1665 


x 10-7Qz? + 77.034LogQz (Range of °C : 25-250) (2) 


Chalcedony 
TCC) = zo Lgo -273.15 (Range of °C : 25-250) 8) 
TC) = 1112 -27315 (Range of °C : 25-250) (4) 
2") = 791-Log Oz 7 8 : 


(b) Na-K-Ca geo-thermometers 
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Na-K 

TCO = oaos togNayR) 273.15 (Range of °C : 25-250) (5) 
Na-C 

Mm1CO= SOR Log ea 273-15 (Range of °C : 25-250) (6) 
Na-K-Ca 

TCC) = 1647. —273.15 (Range of °C : 0-350) 


Log(Na/K)+/Log(Ca?° /Na)+2.24 


(7) 


3.3.2. Ternary diagrams 

Three ternary (CI-SO,-HCO3 and Na+K-Ca—Mg) and Giggen- 
bach (Na-K—Mg) diagrams are constructed for Ain Al Harrah hot 
spring to classify its thermal water (based on major anions and 
cations) and to indicate the prevailing subsurface thermal condi- 
tions [30,31]. 
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3.4. Geothermal reserve estimation 


Geothermal reserves are defined as quantities of thermal 
energy which are expected to be recovered from known reservoirs 
from a given date forward. The reserve is that part of the resources 
which could be extracted economically and legally at present and 
that are known and characterized by drilling or by geochemical, 
geophysical and geological evidence [32,33]. 

The most common approach for geothermal reserve estimation 
is the volumetric method. This method is mainly utilized in the 
later stages of geothermal exploration and after that in the early 
stages of field development to justify drilling and commitment for 
a specified power plant size, although may be used if good surface 
indications and geophysical surveys are available [34]. It can be 
better applied during this stage than numerical modeling which 
requires a significant number of wells and production history to be 
considered reliable [35]. 

In general, the reservoir parameters necessary for geothermal- 
reserve estimation are gathered from the analyses of the tem- 
perature data and from other geophysical and geochemical 
interpretations. A computer program (RESPAR, ICEBOX software) 
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Fig. 9. Accumulation grid map of the intersection, length and frequency of the extratced lineaments of Wadi Al-Lith basin. 
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is used to estimate the different reservoir properties providing 
that the temperature, porosity and the reservoir thickness 
parameters are known. The total stored energy of the geothermal 
reservoir can be regarded as the sum of energy stored in the rock 
(E,) and the energy represented by the fluid stored in the pore 
spaces (Ef). It can be estimated using the following equation: 


E; = Er + Ep = V(1—-#)p,Cr(Ti-To) + Vb pwCw(Ti-To) (8) 


Where £ is the total thermal energy (J) in the rock (E,) and fluid 
(Eş); ¢ is the reservoir porosity (%); V is the reservoir volume (m°); 
Prw are the densities of rock and water (kg/m); C,w are the heat 
capacities of rock and water (J/kg °C) and Ti, T, are the initial 
reservoir and the reference temperatures (°C). 

This stored thermal energy can be converted into power 
potential using the following equation: 


E; x RF x CE 


Power Potential (MWt) = PL LE 
x 


(9) 


where RF is the recovery factor, CE is the conversion efficiency, PL 
is the geothermal plant life in years and LF is the load factor. 
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Fig. 10. Base map of Ain Al Harrah hot spring, showing the location of the 2D 
electric profiles. 
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4. Results and discussion 
4.1, Lineaments extraction and analysis 


The analysis of lineaments has clarified that about 10,000 
elements are detected in Wadi Al-Lith basin (Fig. 4). RockWorks 
software is mainly utilized for drawing the rose diagram for 
detecting the general trends of these lineaments in Wadi Al-Lith 
basin and more specifically in Ain Al Harrah area (Fig. 5). The 
diagram illustrates that the prevailing structural elements are 
nearly the same and almost all lineaments are oriented and 
trended from 0° to 90° (or 180°-270°) with mean trend vector of 
33° (or 213°). Few lineaments are oriented along a northwest- 
southeast direction in Ain Al Harrah hot spring area (Fig. 5b). 
Tables 3 and 4 summarize the main structural trends of these 
lineaments, as well as their percent, lengths and frequency. 

Figs. 6-8 show the constructed intersection, length and fre- 
quency property maps for all the extracted lineaments in Wadi Al- 
Lith area. An accumulation grid map for these three properties is 
drawn, to illustrate their combined effect (Fig. 9). This map shows 
how much these three components behave and affect the study 
area. The weighting factor or influence of each component in the 
accumulation grid map was chosen to be, lineation intersection 
(3), lineation length (2) and lineation frequency (1). 

Interpretation of these figures shows that, no more than three 
lineament elements are acting together in the study area. This 
assigns middle to low influence of these elements in Wadi Al-Lith 
catchment area. Outside Wadi Al-Lith area, on the other hand, a 
high combined influence of these structural lineaments is 
recorded in the order of four elements (Fig. 9). The result that 
can be concluded from this lineament analysis is that the surface 
structural elements are not the factor that control the movement 
of the deep arising thermal water of the main hot spring in Wadi 
Al-Lith area, and that there might be another subsurface elements 
that could carry out this geothermal water to the surface. 


4.2. Interpretation of geophysical profiles 


Out of the nine conducted 2D electric resistivity sections, three 
are selected to demonstrate the prevailing structural elements. 
These profiles are specifically conducted in the area of the hot 
spring (Fig. 10). The results obtained from these profiles are mainly 
concerned with the pathways of geothermal water, detecting feed 
zones, the influencing structural elements and the possible facture 
system. 

Fig. 11 exhibits the interpreted 2D electric section (Line I) that is 
extended along an east-west direction. It shows a tilted thermal 
feed zone at the middle of the section oriented along a very 
obvious fault zone. The strike of this fault is found to be along a 
northwest-southeast direction which is considered the main 
prevailing structure in the area (Red Sea Rift). This structure trend 
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Fig. 11. 2D interpreted resistivity profile (line I) extending East-West, Wadi Al-Lith area. 
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Fig. 12. 2D interpreted resistivity profile (line H) extending East-West, Wadi Al-Lith area. 
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Fig. 13. 2D interpreted resistivity profile (line B) extending North-South, Wadi Al-Lith area. 


Table 5 


Petro-thermal parameters of Ain Al Harrah hot spring, Wadi Al-Lith area. 


Hot spring Subsurface temp. (°C) Discharge enthalpy (kJ/kg) Heat flow (mW/M°) 
Ain Al Harrah 136 219 183 
0 ¢l100 Nat+K 
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Fig. 14. Type of thermal fluids of Ain Al Harrah hot spring as inferred from; (A) CI-SO,-HCO; and (B) NA+K-Ca—Mg ternary diagrams. 


is opposite to another oriented northeast-southwest geothermal 
feed zone, detected at the western end of the profile H between 
electrodes 80-90 (Fig. 12). 


Fig. 13, on the other hand, represents another 2D electric profile 
extending from north to south (Line B). A suggested fault system is 
recognized at the middle of the section. This may be interpreted as a 
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Fig. 15. Giggenbach diagram of Ain Al Harrah hot spring, Wadi Al-Lith area. 


result of the well known northeast-southwest trending faults that 
influence the western coastal parts of the Red Sea area. In general, 
the feed zones of thermal waters are indicated by lower resistivity 
values and the structural control of the thermal water is obvious. 


4.3. Characteristics of thermal water 


The petro-thermal parameters which are indicated from the 
geo-thermometers are listed in Table 5. It appears clearly that the 
Ain Al Harrah hot spring attains good and high geothermal 
parameters. The weighted values of the estimated subsurface 
formation temperature, discharge enthalpy and heat flow are 
found to be 136 °C, 219 kJ/kg and 183 mW/M?, respectively. 

Figs. 14 and 15 represent the constructed ternary (CI-SO,-HCO3 
and Na+K-Ca-Mg) and Giggenbach (Na-K-Mg) diagrams for Ain Al 
Harrah hot spring and one neighboring water well. These dia- 
grams are used to indicate the type of the geothermal water and 
the subsurface geothermal condition at which the dissolved ions 
of the surface ascending thermal fluids are originated. 

The ternary diagrams illustrate that the water type of the hot 
spring is completely different than that of the neighboring water 
wells. The increase of the concentration of the chlorine anion (mg/l) 
on the expanse of the sulfate and carbonate inions is very clear on 
the hot spring (Fig. 14a). On the other hand, a considerable increase 
of the Mg cation is indicated for the water well sample on the 
expanse of the Na+K and Ca cations (Fig. 14b). A mature water type, 
most probably mixed with surface water is indicated to the hot 
spring, while a volcanic unmixed water type is given for the water 
well. Giggenbach (Na-K—Mg) diagram shows that Ain Al Harrah hot 
spring assigns higher geothermal parameters than the water well. 
The cluster of this hot spring is located along the 150 °C Mg-K and 
the 175 °C-225 °C Na-K iso-thermal lines, while that of the water 
well is located very close to the Mg corner along iso-thermal lines of 
less than 50 °C assuming very low geothermal regime (Fig. 15). 

Another unique difference in water type, is the increase of the 
concentration of some radioactive elements in the hot spring like 
rubidium (Rb, 68.9 ppm) and strontium (Sr, 4160 ppm) more than 
those in the water well (Rb, 0.587 ppm and Sr, 4160 ppm). 


4.4. Geothermal energy potential 


In this study, geothermal-reserve estimation is carried out for 
evaluating the possible geothermal capacities of Ain Al Harrah hot 


Table 6 
Reservoir parameters used for geothermal reserve estimation. 


Property Unit Value 
Reservoir Volume km? 1.28 
Porosity % 4 
Density Rock kg/m? 2800 
Water kg/m? 980 
Heat capacity Rock J/kg °C 1000 
Water J/kg °c 4200 
Temperature Initial °C 95 
Reference C 45 


spring. Reserve estimation is based on the above mentioned Eqs. 
(8) and (9), assuming the following parameters: 


e A reservoir area of 4.8 km? (1.2 x 4.0 km?) and reservoir volume 
of 1.28 km? providing a reservoir thickness of 250 m as con- 
cluded from the interpretation of resistivity data. 

è An initial temperature of 95 °C for the geothermal reservoir. 

è A reference temperature of 45°C is suggested (due to the 
higher average annual surface temperature in Saudi Arabia). 

© Pore volume range of 3-5 p.u. for the fractured granitic rocks 
(average value of 4 p.u.). 

© Density range of 2.7-2.90 g/cc (average value of 2.80 g/cc). 


The different reservoir parameters that are derived from the 
analyses of temperature and geo-thermometer data and can be used 
in thermal energy estimation are represented in Table 6. According to 
these parameters and by applying in Eq. (8), the total stored heat 
energy of Ain Al Harrah geothermal area is found to be 1.713 x 10” J 
(1.34 x 10'€ J for thermal fluids and 1.57 x 10” J for reservoir rock). 

The total stored heat energy of Ain Al Harrah hot spring can be 
converted into power energy units (megawatt thermal using Eq. 
(9)) for possible energy production and other low geothermal 
utilizations. Assuming recovery factor of 0.2 (normal factor for 
geothermal reservoirs), turbine conversion efficiency of 0.47, 
(initial temperature/reservoir temperature) plant life of 20 years 
and load factor of 0.95, the thermal power potential of Ain Al 
Harrah geothermal reservoir is found to be 26.99 MWt. This 
geothermal reserve is considered a very promised and can be 
used for providing the Al Lith area with electricity (clean energy) 
on a long term basis. We strongly recommend exploiting this 
geothermal system to generate electricity for the first time in 
Saudi Arabia using a specific power plant that operates thermal 
fluid that boils at low temperature (e.g. Kalina Cycle). 


4.4.1. Utilization of geothermal energy in Saudi Arabia and economic 
and social impact 

The estimated geothermal energy ( MWt) is quite enough for 
possible electricity production in Saudi Arabia. Kalina power plants 
can be used due to the lower boiling point (70 °C) of their fluids. 
Having a higher temperature range than the boiling point of the 
Kalina cycle, the thermal water of Ain Al Harrah hot springs can be 
economically used for power generation. However, this thermal 
water can be utilized directly in other low forms of applications 
(district heating, fish farming, agricultural applications, greenhouses, 
etc.). Some of these applications are already now in use (Gouba area, 
central Saudi Arabia), where a number of geothermal-based natural 
therapy, medical and swimming pools are constructed. 

The social factor of using the geothermal energy implies from 
the fact that most of the geothermal resources are mainly located 
in the extreme western and south-western borders of Saudi Arabia 
along the coast of the Red Sea. These areas besides being far are 
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topographically complicated, not so easy to reach and need more 
developments in its infrastructures and the level of the residence 
Bedouins people. The geothermal based-energy can help in con- 
structing new modern communities fed with the geothermal- 
based energy, where no high voltage transmission lines are 
available in the vicinity and where it would be too expensive to 
connect to the national electric grid. 


5. Conclusions 


The following conclusions may be drawn according to the 
present study: 


1. Wadi Al-Lith area contains one of the most promised geother- 
mal systems in the western of Saudi Arabia. It is located near 
the coast of the Red Sea and is strongly affected by its tectonic 
activity and structural regimes. 

2. Ain Al Harrah hot spring is considered the main geothermal 
target in Wadi Al-Lith area that attains high surface tempera- 
ture upto 95°C and good petro-thermal characteristics of 
136 °C (reservoir temperature), 219 kJ/kg (discharge enthalpy) 
and 183 mW/M? (heat flow) as indicated from the geo- 
thermometric analysis. 

3. The analysis of Landsat RGB images (lineaments extraction) and 
the interpretation of the electric geophysical profiles have 
clarified the presence of many surface and subsurface struc- 
tural elements. Feed zones of thermal waters are indicated by 
lower resistivity values and the structural control of the 
thermal water is obvious. Majority of the encountered 
geothermal feed zones are tilted and aligned along a north- 
east-southwest (transform) direction and little bit to a north- 
west-southeast one (Red Sea Rift). 

4. Almost all surface structural lineaments that control the move- 
ment of the geothermal water (Ain Al Harrah hot spring) are 
found following the main structural trends of Wadi Al- 
Lith area. 

5. The analysis of ternary and Giggenbach diagrams shows that 
the water type of geothermal system is completely different 
than the neighboring domestic water aquifers. A mature mixed 
water type is assigned for Ain Al Harrah hot, while a volcanic 
unmixed water model is given for the neighboring water wells. 

6. The stored heat energy is found to be 1.34 x 10'°J for the 
thermal fluids and 1.57 x 10!” J for the hot rock. 

7. A geothermal power potential of 26.99 MWt is estimated for 
Ain Al Harrah hot spring which is considered good and 
adequate for economic investment. 

8. A medium scale power plant that utilizes fluids of low boiling 
point is recommended for possible electricity generation for 
the first time in Saudi Arabia. 

9. A very strong and positive economic-social impact will be 
obtained upon applying such geothermal technology in the 
study area or in other remote areas of similar geothermal 
systems. It will enhance the infrastructure and build up a 
new geothermal-feed communities based on clean renewable 
energy. 
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